Exact eigenvalue spectrum of a class of fractal scale-free networks 
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Abstract - The eigenvalue spectrum of the transition matrix of a network encodes important 
information about its structural and dynamical properties. We study the transition matrix of 
a family of fractal scale-free networks and analytically determine all the eigenvalues and their 
degeneracies. We then use these eigenvalues to evaluate the closed-form solution to the eigentime 
for random walks on the networks under consideration. Through the connection between the 
spectrum of transition matrix and the number of spanning trees, we corroborate the obtained 
eigenvalues and their multiplicities. 
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Introduction. — Complex networks have emerged as 
a universal tool for studying complex systems in many 
different fields, such as physics, chemistry, biology, and 
computer science [TH3]- A fundamental issue in network 
science is to understand diverse dynamical processes de- 
fined on different networks |3HS]. A paradigmatic dynam- 
ical process that has attracted increasing research interest 
is random walks [?Hl 7j. which have found a broad range 
of applications in various areas of science and engineer- 
ing HHH2H- 

As is well known, random walks on a network can be 
described by a transition matrix and various interesting 
quantities about random walks are closely related to the 
spectrum of the transition matrix [IT1I22H24] . For exam- 
ple, the mean first-passage time (MFPT) from one node 
to another on a network can be expressed in terms of the 
eigenvalues and orthonormalizcd eigenvectors of its tran- 
sition matrix [XT] [221 [23]; the sum of reciprocals of one 
minus each eigenvalue, other than the eigenvalue f, de- 
termines the eigentime identity [23| , sometime referred to 
as Kemeny's constant [21]; the second largest eigenvalue, 
together with the smallest eigenvalue, determines the mix 
time of random walks on the network |22[|23| . Thus, it is 
of theoretical and practical significance to calculate and 
analyze the spectra of the transition matrices of networks. 
However, little attention has been paid to this relevant is- 
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sue, the main reason for which is that the computational 
complexity for computing the eigenvalue spectrum of a 
general network is very high |25j . 

In this paper, we present a theoretical study of the tran- 
sition matrix for unbiased random walks on a family of 
iterative networks [2l)rf2"5|. which have remarkable scale- 
free [22] and fractal [3U] phenomena that constitute our 
fundamental understanding of real-life systems. Making 
use of the real-space decimation technique [3T1I32] , we will 
derive an exact recursion expression for the eigenvalues at 
every two consecutive iterations, which is obtained based 
on the structure of the studied networks. We then proceed 
to provide the complete spectrum of eigenvalues and their 
corresponding degeneracies. Using these eigenvalues, we 
further derive an explicit formula for the eigentime iden- 
tity. Finally, we prove the validity of the computation for 
eigenvalues based on the relationship between the num- 
ber of spanning trees and the eigenvlaues of the transition 
matrix of a network. 

Model and properties. — The considered frac- 
tal scale-free networks arc constructed in an iterative 
way [261428J . Let F n (n > 0) denote the networks after 
n iterations. For n — 0, Fq consists of an edge connecting 
two nodes. For n > 1, F n is obtained from F n -\ by per- 
forming the following operations on every existing edge in 
F n -\ as shown in Fig. Q] replace the edge by a path of 2 
links long, with both endpoints of the path being the iden- 
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Fig. 1: (Color online) Construction of the networks. One 
can obtain the next generation of the network family through 
replacing each edge of the present generation by the clusters 
on the right-hand side of the arrow. 
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tical cndpoints of the original edge, then for each endpoint 
of the path create m (a positive integer) new nodes and 
attach them to the endpoint. Figure [5] illustrates the first 
several iterative construction steps for a particular case 
of m = 1, whose dynamics of random walks were studied 
in J33]. 

The deterministic construction of the networks allows 
for treating exactly their relevant properties. At each gen- 
eration rii (rii > 1), the number of newly introduced nodes 
is V ni = (2m + l)(2m + 2)" i ~ 1 . Then, the network size 
(number of nodes) N n of F n is given by 



N„ 



(2m + 2) 



1. 



(1) 



And the total number of edges is 

E n =V n -l= (2m + 2)™ . (2) 

Let di(n) be the degree of node i in F n that entered the 
networks at generation rii (n-i > 0). Then di(n + 1) = 
(m + 1) di(n), that is, after each new iteration the degree 
of node i increases by m times. 

These resultant networks display the following inter- 
esting structural features. They are in power law with 
the degree distribution exponent being 7 = 1 + ln(2?n + 
2) / ln(m + 1) [27ll34U35j . In addition, after the evolution of 
one generation, the number of nodes increases by a factor 
Jn = 2m + 2, see Eq. (fTJ; the diameter grows by a factor 
of fo = 2 [57] , and the MFPT between two previously ex- 
isting nodes increases by a factor of /m = 4m + 4 [2"TIl3"5] . 
Thus, the networks are fractal with a fractal dimension 
fs = ln/jv/ln/jj = ln(2m + 2)/ln2 [HI [26l [34] 5 their 
random- walk dimension is f w = In Jm / In /d = ln(4m + 
4)/ In 2, and their spectral dimension is f s = 2fs/f w = 
2 ln(2m + 2)/ ln(4m + 4). Finally, they are "large- world" 
with their diameter and average distance growing in a 
power of the network size [27ll34l[35] . 



Spectrum of transition matrix. — After introduc- 
ing the model and some properties of the networks, in this 
section we study the eigenvalue spectrum of the transition 
matrix for the network family, as well as their degenera- 
cies. 

Eigenvalue spectrum. As is well known, the structure 
of F n is encoded in its adjacency matrix A n , whose en- 
tries A n (i,j) are defined by A n (i,j) = 1 if nodes i and 



Fig. 2: (Color online) Growth process for a special case of 
network corresponding to m = 1. 



j are adjacent in F n , or A n (i,j) = otherwise. Then, 
the transition matrix of F n , denoted by T„, is defined as 
T„ = D^An, where D n is the diagonal degree matrix of 
F n with its iih diagonal entry being di(n). Thus, the el- 
ement of T n is T n (i,j) = A n (i, j)/di(n), which represents 
the jumping probability of the discrete-time unbiased ran- 
dom walks [7J for a particle going from node i to node j. 

We now consider the eigenvalue spectrum of T n . Since 
T n is asymmetric, we introduce the following matrix: 



DnT n D It 



(3) 



Obviously, P n is real and symmetric and has the same set 
of eigenvalues as T n . It is easy to verify that the entry of 

= . Next, we apply the decima- 



P n is P n (i,j) 

tion method 31 , 32] to determine the eigenvalues and their 
multiplicities of P n . The decimation approach is universal 
and has been used to compute the Laplacian spectra of 
Vicsek fractals [36H38] and their extensions [39] . 

We now address the eigenvalue problem for matrix 
P n+ \. Let a represent the set of nodes belonging to F n , 
and (3 the set of nodes generated at (n + l)th iteration. 
By definition, P„+i has a block form 



Pn 



+ 1 



P n 



P 



0,a 



Pf>,P 











(4) 



where we have used the facts that P a , a is the zero matrix 
of order N n x N n since there are no transitions between 
any pair of original nodes, and that Pp,p is the zero matrix 
of order (N n+1 - N n ) x (N n+1 - N n ). 

Suppose Aj(n + 1) is an eigenvalue of P n +i, and u = 
(u a ,up) T is its associated eigenvector, where the super- 
script T represents transpose and u a and up correspond 
to nodes in a and /?, respectively. Then, one can write 
the eigenvalue equation for matrix P„+i in the following 
block form: 



P a ,p 

P ,a 



Up 



Xi(n + 1) 

Equation ([5]) can be recast into two equations 



up 



P a ,pup = \ (n+l)u a , 



(5) 



(6) 
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P/3, a u a = Xi(n + l)u . 
Equation (J7J yields 

1 



up 



Xi(n + 1) 



(T) 



(8) 



provided that A,(n + 1) ^ 0. Plugging Eq. flgj) into Eq. © 
gives 

Y7 — —rrP a ,pPp, a u a = X,{n + l)u Q . (9) 
Aj + 1 ) 

Thus, we reduce the problem of determining the eigenval- 
ues for matrix P n +i of order N n+ i x N n+ i to calculating 
the eigenvalues of matrix P a ,pPfj. a with a smaller order of 
N n x N n . 

In the Appendix, we prove that 



2m + 1 



1 



2m + 2 



p 



(10) 



where /„ is the identity matrix of the same order as P n . 
Thus, we have related P a ,pPp, a to P n , which enables us to 
express the eigenvalues of matrix P n +i in terms of those 
ofP n . 

Substituting Eq. (fTUft into Eq. ©, we obtain 



1 



2m + 1 



In 



1 



P„ M a = Ai(n + l)u a , 



(11) 



Ai(n + 1) \2m + 2 " ' 2m + 2" 
that is, 

P n u a = {(2m + 2)[Xi(n+ l)] 2 - (2m + 1)} u a . (12) 



Hence, if Aj(n) is the eigenvalue of P ra corresponding to 
the eigenvector u a , Eq. (fT2j) implies 



Ai(n) = (2m + 2)[A*(n + l)] 2 - (2m + 1). (13) 

Solving the quadratic equation in the variable Xi(n + 1) 
given by Eq. (fT3|) yields 



A^n+l) = V ^^V" ' A ^ 2 ^ + 1 ) = - A m(" + !) • 

(14) 

Equation (fl"4"|) relates Aj(n + 1) to Aj(n), with each Aj(n) 
giving rise to two eigenvalues of P„+i. 

Degeneracies of eigenvalues. To determine the multi- 
plicities of the eigenvalues, we first calculate numerically 
the eigenvalues for those networks of small sizes. For Fq, 
the eigenvalues are 1 and — 1; while for Pi, its eigenval- 



ues are 1, — 1, 0, 



and 



-Vm+7; 



For n > 2, we 



m+l ' m+1 

find that the eigenvalue spectrum displays the following 
properties: (i) From Eq. (|T4]) . eigenvalue 1 gives rise to 
two eigenvalues 1 and —1 with single degeneracy, (ii) All 
eigenvalues of a given generation n% always exist at its 
subsequent generation n, + 1, and all new eigenvalues at 
generation n» + 1 are just those generated via Eq. (fl"4"|) 
by substituting Xi(n) with Xi(n{) that are newly added 
to generation nf, moreover each new eigenvalue keeps the 



degeneracy of its father, (hi) Except for 0, all other eigen- 
values are generated from —1 and 0, and all the offspring 
eigenvalues of —1 are nondegenerate. Thus, all that is left 
is to determine the degeneracy of and the multiplicities 
of its offsprings, based on property (ii). 

Let .D™ ul (A) denote the multiplicity of eigenvalue A of 
matrix P n . We now find the number of eigenvalue of P n . 
Let r(M) be the rank of matrix M. Then, the degeneracy 
of eigenvalue of P„+i is 



D™[(X = 0) = N n+1 -r(P n+1 ). 



(15) 



In order to determine D^Y\(X = 0), we can alterna- 
tively compute r(P„+i). Obviously, r(P n +i) = r(P a ,p) + 
r{Pfi, a ) = 2r(P (3 ^), where r(P a> p) = r(Pp >a ) is used. 

We next determine r(P^ iQ ). First, we show that Pp, a is 
a full column rank matrix. Let 



N„ 



-nJ T = Y,kiMi, (16) 



where Mi is the column vector of Pp, a as defined in 
Eq. (061). Let Mi = {M hl ,M 2il ,..'.,M Nn+1 - Nnil ) T . 
Suppose that v = 0. Then, we can prove that for 
an arbitrary hi, ki = always holds. By construc- 
tion, for any old node i £ a, there exists a new leaf 
node / <G (3 attached to i. Then, for vi = kiM\ t i + 
hM 2 ,i, • • • , k Nn+1 -N n ,iM Nn+1 -N n ,l, only M M ^ but all 
M x ,i = for x 7^ i. From vi = 0, we have ki = 0. There- 
fore, r(Pp tCC ) = N n = (2m + 2)" + 1. 

Combining the above-obtained results, we obtain the 
degeneracy of eigenvalue of P n as 



L>r'(A = 0) 



0, 

2m(2m + 2) ? 



n = 0, 
1, n>l. 



(17) 



Since every eigenvalue in P n keeps the degeneracy of its fa- 
ther, the multiplicity of each first-generation descendants 
of eigenvalue is 2m(2m + 2)"~ 2 — 1, the multiplicity 
of each second-generation descendants of eigenvalue is 
2m(2m + 2)"~ 3 — 1, and so on. Thus, the total number of 
eigenvalue and all of its descendants in P„ (n > 1) is 



A^ ccd (0) = ^[2m(2m + 2)'*- 1 - 1]2" 



(2 m + 2 



|7i ry7l-\-l 



1. 



(18) 



Similarly, the total number of eigenvalue —1 and its de- 
scendants in P„ (n > 0) is 



K ccd (-i) = 21 = 2 " +1 ~ L 



(19) 



i=0 



Summing up the number of eigenvalues found above, we 
have 

A^ ocd (0) + A^ cod (-l) + 1 = (2m + 2)" + 1 = N ni (20) 

which means that we have determined all the eigenvalues 
of P„. 
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Application of eigenvalue spectrum. We next 
show how to use the obtained eigenvalues and their multi- 
plicities to determine some related quantities for the frac- 
tal scale- free network family F ni including the eigentime 
identity of random walks and the number of spanning 
trees. We note that the number of spanning trees for every 
generation is 1, so our aim for evaluating spanning trees 
is to verify that our computation of the eigenvalues of the 
transition matrix is right. 

Eigentime identity. Eigentime identity of random 
walks on a network is a global characteristic of the net- 
work, which reflects the structure of the whole network. 

Let Hij(n) denote the mean-first passage time from 
node i to node j in F n , which is the expected time for 
a walker starting from node i to reach node j for the first 
time. Let it = (tti, tt2, . . . , 7r/v) T represent the stationary 
distribution for random walks on F„ [22|[23] . It is easy 



to derive iTi — di(n)/(2E n ), satisfying Y^iLi^i = 1 an d 
7r T T„ = 7r T . Then, the eigentime identity, denoted by 
H n , for random walks on F n , is the expected time from 
a node i to another node j, selected randomly from all 
nodes accordingly to the stationary distribution. That is, 



N„ 



1 



(21) 



where we have assumed that Xi(n) = 1. 

Note that for any i, Ui{n) = l — \i(n) is an eigenvalue of 
the normalized Laplacian matrix of F n , defined as L n = 
In — Pn [40H42] . The normalized Laplacian matrix has 
found numerous applications [431445] . By definition, 



Ft,, 



E 



i 



<Ti(n) ' 



(22) 



The goal next is to explicitly evaluate this sum. 

Let fi n be the set of all the N n — 1 nonzero eigenvalues 
of matrix L„, fi„ = {02(71), ^s{n), &N n {n)}, in which 
the distinctness of the elements has been ignored. It is 
clear that fi n (n > 1) includes 1, 2, and other eigenval- 
ues generated by them. Then, fi„ can be classified into 
three subsets represented by fi n , fin and fi n , respec- 



tively. That is, fi„ = fi^ 1 U Vl)? 1 U fin J; , where fi^ con 
sists of eigenvalue 1 with multiplicity 2m(2m + 2) r 



(2) 



,(3) 



,(1) 



-1, 



(2) 

fin contains only eigenvalue 2 with a single degeneracy, 

(31 

and fi„ includes those eigenvalues generated by 1 and 



2. Obviously, ^\ 

E 



ie£)Jf J <Ji(n) 



ieiV„> <7i(n) 



2m(2m + 2)' 



1 and 



Each eigenvalue Ui{n — 1) in fi^ 3 ^ generates two eigen- 

values, CTi ; i(n) and <Ji^{n), belonging to fin , via the fol- 
lowing equation: 



[CTi(n)] 2 - 2cr. t (n) + 



Qj{n - 1) 
2(m + 1) 



0. 



(23) 



which is easily obtained from Eq. (fT3"]l. According to 
Vieta's formulas, we have <Ji,\(n) + <Jifl.(j>) = 2 and 
(Jis(n) x cr i;2 (n) = 2(m+i) ■ Tnen 



1 



1 



4(m+ 1) 



<?*,i(n) o'i,i( n ) <7j(n— 1)' 
which implies that 

1 _ NT- 1 



E 



o-i(n) 



4(m + 1) 



Thus, we have 



<Ti(n-i)en 



(3) 



o-i(n - 1) 



H n = 4(m + l)ff„_i + [2m(2m + 2)' 



- 11 



(24) 



(25) 



(26) 



Considering the initial condition Hi = 8rr ]+ 3 1 Eq. ([26]) can 
be solved by induction to yield 

(6m 2 + 6m + l)(4m + 4)" 



H,, 



(m + l)(4m + 3) 



-2m(2m + 2 



n-l 



1 



8m + 6 

which can be expressed in terms of network size N n as 

6m 2 + 6m + 1 _ 1 %l+ln2/ln(2m+2) 

(m + l)(4m + 3r J 



(27) 



■(JV n 



1) 



1 



m + 1 8m + 6 

Thus, for large networks, i.e., A 9 — > 00, 

H n ~ (iY ri )l+'n2/ln( 2m+ 2) = (^^I + V/b 



(28) 



(29) 



increasing as a power-law function of the network size N n 
with the power exponent greater than 1. Note that the 
power-law exponent in Eq. ([2^]) is consistent with the gen- 
eral scaling for the average global MFPT (i.e., eigentime 
identity) previously obtained in [13] , but is comparable to 
the diverse behaviors of MFPT observed in real-life sys- 
tems, the reason for which was explained in |46j . 

Spanning trees. An important task in the study of 
networks is to determine the number of spanning trees on 
different networks [47HBT] . A spanning tree of a connected 
graph is defined as a maximal set of edges of the graph that 
contains no cycle, or as a minimal set of edges that connect 
all nodes. It has been shown [40l|42] that for a connected 
network G with N nodes, the number of its spanning trees, 
N s t(G), is related to the N — 1 nonzero eigenvalues of its 
normalized Laplacian matrix. According to [101113, the 
number of spanning trees N s t(F n ) for F n is determined by 



N 1 1 



N 8t (F n ) = [piwn^w / E d *( n )- ( 3 °) 

i=l i=2 I i=l 

We now determine the terms in the sum and in the 
products of Eq. ([30]) . It is evident that 



di(n) = 2E n = (2m + 2)". 



(31) 
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Let A n and A„ denote Yli=i di{n) and Yli^2 a i( n )> respec- 
tively. According to the above-obtained results, it is easy 
to show that they obey the following recursive relations: 



A n _i x (m + 1)^—1 x 2 



and 



A, 



1 



2m + 2 



A, 



(32) 



(33) 



Using Ai = 1 and Ai = 2, we can solve Eqs. (j3"2"|) and (|33l) 
to obtain 



A„ = (to + l)"(2m + 2) <2 "™+i" 



and 



A„ 



_ nn+l 



1 



2m + 2 



(2m + 2)" -1 
2m + l 



(34) 



(35) 



Inserting Eqs. ([31]), ([M]) and ([35]) into Eq. ([201) gives 
N s t(F n ) = 1, which proves that the technique and process 
of our computation on the eigenvalues and their degenera- 
cies for the transition matrix of the network family F n are 
indeed correct. 

Conclusions. — Rich information about the structure 
and random- walk dynamics of a network can be extracted 
from the spectrum of the transition matrix of the network. 
In this paper, we have studied the transition matrix of a 
class of networks with scale- free and fractal behaviors that 
are observed in various real- world networked systems. The 
iterative construction of the networks allows a detailed 
analysis of the eigenvalues and their multiplicities. Using 
the renormalization approach, we have derived an explicit 
recursive relationship between the eigenvalues of the net- 
work family at two successive iterations. For each eigen- 
value, we have also determined its degeneracy. Based on 
the obtained recursion relation of eigenvalues, we have fur- 
ther determined the eigentime identity for random walks, 
an important quantity rooted in the inherent topology of 
a network. Moreover, we have tested our computation 
and results for the spectrum of the transition matrix via 
evaluating the number of spanning trees in the network 
family under consideration, using the established relation 
between the two quantities. Our work opens up the pos- 
sibilities to characterize the spectrum for the transition 
matrices of other deterministic networks. 

It deserves to be mentioned that by definition eigentime 
identity is the expected time from a node i to another node 
j that is chosen randomly from all nodes according to the 
stationary distribution, which means that the target is im- 
mobile. Then, interesting questions arise: What happens 
if the target is not immobile? Would it be possible in 
this case to use the so-called Pascal Principle [13,52,53 
to obtain a rigorous upper bound? Although the ques- 
tions go beyond the scope of the present paper, they are 
worth studying in the future. Futhermore, future work 
should also include how to apply the theory to real- world 
networks, a relation between whose topology and function 
has been demonstrated l54l. 
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Proof of equation ()10p . — In order to prove 
the equivalent relation P a .pPp^ a = j^I n + 2^+2 P n> 
it suffices to show that the entries of the former are 
equal to their counterparts of the latter. For matrix 
2m+2 ^ n + 2m+2 Pn ' denoted by Q„, obviously its entries 
are: Q„{i, i) = §^±± for i ^ j and Q n (i,j) = 5^+2-^ (*>i) 
otherwise. 

We proceed to determine the entries of P a ,pPp,on de- 
noted by R n . Note that P a ,p can be written as 



( ^ \ 



P. 



Mj 



(36) 



where each Mj is a column vector with order N n+ i — N n . 
Since P aJ3 = Pj a , we have P^ a = (Mi M 2 ...M Nn ). 
Then, the entry R n (i,j) = Mj Mj of R n can be deter- 
mined in the following way. 

If % = j, the diagonal element of R n is 



Rn(i,i) = MjM, 

^2 P n+1 (i, k)P n+1 (k,i) = ^2 

fce/3 fce/3 



A n+ i(i, k) 



di(n + l)d k (n + 1) 



1 



di(n + 1) 
1 



E 1+ E 



\ 



di(n+l) 
2m + 1 

2mT^ =Q " M 



fce/3,i~fc fcG/3,i~fc 
\d k (n+l)=2 d k (n+l) = l ) 



di(n) x - + mdi(n) 



(37) 



where the relation d,(n + l) = (m + 1) di(n) has been used. 
In Eq. (|37[) . i ~ k means that nodes i and k are adjacent 
in F n+X . 

If i 7^ j, the non-diagonal element of R n is 

R n {i,j) = MjM j 
= ^2Pn+i(i,k)P n +i(k 1 j) 



E 



A n+1 (i,k) 
A„ +1 (fe,j)=l 



1 d k (n+ l) v / d l (n + l)d j {n + l) 
1 



(2m + l)y/ di(n)dj(n) 
Thus, Eq. (|TU1) is proved. 



Qn(iJ)- 



(38) 
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